Sharp increase of the effective mass near the critical density 
in a metalhc 2D electron system 
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We find that at intermediate temperatures, the metaUic temperature dependence of the conduc- 
tivity (t(T) of 2D electrons in siUcon is described well by a recent interaction-based theory of Zala 
et al. (Phys. Rev. B 64, 214204 (2001)). The tendency of the slope a~^da/dT to diverge near 
the critical electron density is in agreement with the previously suggested ferromagnetic instability 
in this electron system. Unexpectedly, it is found to originate from the sharp enhancement of the 
effective mass, while the effective Lande g factor remains nearly constant and close to its value in 
bulk silicon. 
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In dilute two-dimensional (2D) electron systems, the 
energy of electron-electron interactions dominates the ki- 
netic energy making the system strongly correlated [Q. 
The interaction strength is characterized by the Wigner- 
Seitz radius, r^, which is equal in the single- valley case to 
the ratio of the Coulomb and the Fermi energies. Accord- 
ing to the Fermi liquid theory |^] (whose applicability to 
dilute 2D electron systems is discussed in Ref. [||), the 
electron-electron interactions should give rise to a renor- 
malization of the system parameters including the effec- 
tive electron mass, m, and the effective g factor. A sharp 
enhancement of the product gm with decreasing electron 
density — possibly, a precursor of the ferromagnetic in- 
stability — has been observed in recent studies of the par- 
allel field magnetoresistance of a metallic 2D electron sys- 
tem in high-mobility silicon metal-oxide-semiconductor 
field-effect transistors (MOSFETs) [§. This agrees well 
with the gm value obtained from the beating pattern of 
Shubnikov-de Haas oscillations in tilted magnetic fields 

Hi- 

Recently, temperature-dependent corrections to con- 
ductivity due to electron-electron interactions have been 
calculated by Zala et al. ^ based on the Fermi liquid 
approach. In contrast to pre-existing theories (see, e.g., 
Refs. 1^,0), the new theory incorporates strongly inter- 
acting 2D electron systems with electron densities, Us, 
down to the vicinity of the metal-insulator transition 
(provided that the conductivity a ^ e^/h). For suffi- 
ciently strong interactions, it predicts a metallic temper- 
ature dependence of conductivity in the entire tempera- 
ture range. At very low temperatures, in the "diffusive" 
regime (T ^ h/kgT, where r is the elastic relaxation 
time), this is Finkelstein's weakly- metallic (logarithmic) 
conductivity p^ , p^ . At intermediate temperatures, in 
the "ballistic" regime (T > Ti/fesr; T > 0.2 - 0.5 K 



under the conditions of our experiments), the predicted 
cr(r) is similar to the Gold-Dolgopolov linear dependence 



1 - AkeT, 
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where the slope. A, is determined by the interaction- 
related parameters: the Fermi liquid constants, Fq and 
Ff. These parameters are responsible for the renormal- 
ization of the g factor and the effective mass 0] 
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and can be determined experimentally |l4[ . The slope A 
is predicted to rise as the interaction strength increases 
and the 2D electron system is driven toward the ferro- 
magnetic instability. The last is expected to occur, in 
the simplest case, at Fg = — 1, which corresponds to the 
diverging effective g factor. 

In this paper, we perform precision measurements of 
the temperature-dependent conductivity in a metallic 2D 
electron system in silicon over a wide range of electron 
densities above the critical electron density, ric, for the 
metal-insulator transition. The theory of Zala et al. |^ 
is found to be successful in interpreting the experimental 
data in the ballistic regime. Knowing the product gm 
from independent measurements, we determine both g 
and TO as a function of rig from the slope of the temper- 
ature dependence of the conductivity. The tendency of 
the slope to diverge near the critical density is consis- 
tent with the suggested ferromagnetic instability in this 
electron system P,p^. However, unlike in the simplest 
scenario for the ferromagnetic instability, it is the value of 
the effective mass that becomes strongly enhanced with 
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FIG. 1. The temperature dependence of the normahzed conductiv- 
ity at different electron densities (indicated in units of 10^^ cm~^) 
above the critical electron density for the metal-insulator transi- 
tion. The dashed lines are fits of the linear interval of the depen- 
dence. 

decreasing electron density, while the g factor remains 
nearly constant, g ~ = 2 in bulk silicon. 

Measurements were made in an Oxford dilution refrig- 
erator with a base temperature of « 30 mK on high- 
mobility (lOO)-silicon samples similar to those previously 
used in Ref . . The resistance was measured by a stan- 
dard 4-terminal low-frequency technique. Excitation cur- 
rent was kept low enough to ensure that measurements 
were taken in the linear regime of response. Contact 
resistances in our samples were minimized by using a 
split-gate technique that allows one to maintain a high 
electron density in the vicinity of the contacts regardless 
of its value in the main part of the sample. In this pa- 
per we show results obtained on a sample with a peak 
mobility close to 3 m^/Vs at 0.1 K. 

Typical dependences of the normalized conductivity on 
temperature, a{T)/ao, are displayed in Fig. |l|at different 
electron densities above the critical electron density for 
the metal-insulator transition which in this sample oc- 
curs at Uc — 8 X 10^° cm^^ |l7|; the value ctq, which has 
been used to normalize ct, was obtained by extrapolat- 
ing the linear interval of the (t(T) dependence to T = 0. 
As long as the deviation |(t/(To — 1| is sufSciently small, 
the conductivity a increases linearly with decreasing T 
in agreement with Eq. until it saturates at the lowest 
temperatures |T8| . 

In Fig. ^, we show the Ug dependence of the inverse 
slope 1/A extracted from the (j{T) data. Also shown for 
comparison is the magnetic energy, hbBc = nh^ns/gm 
(where /is is the Bohr magneton), corresponding to the 
onset of full spin polarization in this electron system, 
which is governed by the (enhanced) product gm [Q. 
Over a wide range of electron densities, the values 1/A 
and hbBc turn out to be close to each other. The low 



FIG. 2. Comparison of the inverse slope 1/A (dots) and the data 
for the polarization field B(, (diamonds) as a function of electron 
density. The dashed lines are linear fits which extrapolate to the 
critical electron density for the metal-insulator transition. 

density data for 1/A are approximated well by a linear 
dependence which extrapolates to the critical electron 
density Uc in a similar way to the behavior of the po- 
larization field Be (note that the tendency of the po- 
larization field to vanish near ric was also reported in 
Ref. ||l^). We emphasize that it has been verified with 
the help of the weak-field low-temperature Hall effect 
measurements that the density of the delocalized elec- 
trons in the metallic phase is practically coincident with 

As has already been mentioned, the coefficient A in 
the linear-in-T correction to conductivity of Eq. (|l|) is 
determined by the Fermi liquid parameters 

^^Jl + aFS)9m_ 
nh Us 

The factor a is equal to 8 in our case where the temper- 
ature is small compared to the valley splitting ||l9| . This 
theoretical relation allows us to determine the many- 
body enhanced g factor and mass m separately using the 
data for the slope A and the product gm as a function of 
rig (the latter is known from independent measurements 
similar to those described in Ref. 0|). 

In Fig. 1^, we show the so-determined values g/ga and 
m/rrib as a function of the electron density (the band 
mass mij is equal to 0.19r7ie where me is the free electron 
mass). Note that in the range of Ug studied here, the 
low-temperature conductivity a > Se^/h. The behavior 
of g and m at electron densities below = 3 x 10^^ cm~^ 
(corresponding to « 4.8) turns out to be very different 
from that at electron densities above this value. In the 
high Ug region (lower r^), the enhancement of both g 
and m is relatively small, both values slightly increasing 
with decreasing electron density in agreement with earlier 
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FIG. 3. Renormalization of the effective mass (filled squares) and 
g factor (dots) as a function of electron density. The data from 
Ref. are shown by open symbols. The dashed lines are guides 
to the eye. The inset compares the theoretical dependence of the 
renormalization parameter Fq on Ug (solid line) with the data 
(dots) calculated using Eq. (^) from our g values. 

data [Q. Also, the renormalization of the g factor is 
dominant compared to that of the effective mass, which 
is consistent with theoretical studies [Q. The depen- 
dence g{ns) is described reasonably well by the theory: 
the inset of Fig. ^ compares the theoretical renormal- 
ization parameter Fq = —rs/2{2rs + \/2) ^ to that 
calculated using Eq. (^) and the data for g{ns). 

In contrast, the renormalization in the low Ug (criti- 
cal) region, where Vs 3> 1, is much more striking. As 
the electron density is decreased. Fig. ^ shows that the 
renormalization of the effective mass overshoots abruptly 
while that of the g factor remains relatively small, g ~ go, 
without tending to increase. Hence, the current analysis 
indicates that it is the effective mass that is responsible 
for the drastically enhanced gm value near the metal- 
insulator transition. 

The present results for the effective mass and g fac- 
tor in the critical region can be compared to the data of 
Ref. jsj obtained by analysis of the Shubnikov-de Haas 
oscillations in high-mobility Si MOSFETs. As seen from 
Fig. ^, data obtained from these two different methods 
are similar. Thus, the Fermi-liquid-based theory |^ is ad- 
equate in describing the properties of dilute 2D electron 
systems. 

It is important to discuss another consequence of the 
theory Q: the slope A of the temperature dependence 
of the conductivity should increase as the ferromagnetic 
instability in a dilute 2D electron system is approached. 
Since renormalization parameters have not been theo- 
retically calculated in the limit Ts ^ 1, the simplest 
scenario of the ferromagnetic instability is prompted by 
Eq. (|^): Fq —1 causes the effective g factor (and the 



FIG. 4. The elastic relaxation time versus electron density at a 
temperature of 0.1 K, assuming m = mi, (dots), and taking into 
account the renormalization of m (squares). The dashed lines are 
guides to the eye. 

slope A) to diverge. Experimentally, the slope A tends to 
diverge near the critical electron density in a way similar 
to the behavior of the product gm seen in Fig. ^. This 
is consistent with the conclusion of Ref. Q about the 
possibility of ferromagnetic instability in this electron 
system. At the same time, the simplest scenario of a di- 
verging g factor is not the case; instead, it is the growing 
effective mass which controls the anomalous behavior of 
the dilute 2D electron system near the metal-insulator 
transition. 

The effective mass enhancement was traditionally con- 
sidered to be small, and, therefore, the value m ~ mj, 
was used to calculate some of the important system pa- 
rameters, e.g., the elastic relaxation time t extracted 
from mobility. In Fig. ^, we compare the so-determined 
T (circles) with that calculated taking into account the 
enhancement of the effective mass (squares). As seen 
from the figure, in the range of electron densities stud- 
ied, the corrected t keeps increasing down to the low- 
est rij. Therefore, the mobility drop at low Ug in high- 
mobility Si MOSFETs (see, e.g., Ref. (H) turns out to 
come from the m enhancement rather than from the de- 
crease in T, although the value r is still expected to van- 
ish in the insulating phase. The observed behavior of r 
is consistent with that of the temperature range corre- 
sponding to the ballistic regime (see Fig. |l]), which gives 
additional confidence in our analysis down to the vicin- 
ity of the metal-insulator transition. Finally, values of 
T much larger than those previously estimated yield ap- 
preciably smaller quantum level widths in perpendicu- 
lar magnetic fields, which helps to understand why the 
Shubnikov-de Haas oscillations survive near the metal- 
insulator transition, as well as the origin of the oscilla- 
tions of the metal-insulator phase boundary as a function 
of (perpendicular) magnetic field p3]. 
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In summary, we have studied the temperature- 
dep endent conductivity in a wide range of electron 
den sities above the critical electron densitv for the 
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metal-insulator transition. Using the recent theory of 
interaction-driven corrections to conductivity we 
extract Fermi-liquid parameters from the experimental 
data and determine the many-body enhanced g factor 
and the effective mass. The tendency of the slope A of 
the temperature dependence of the conductivity to di- 
verge near the critical density is in agreement with the 
suggested ferromagnetic instability in this electron sys- 
tem Q. Unexpectedly, it is found to originate from the 
growing effective mass rather than the g factor. In addi- 
tion, the mass enhancement is found to be responsible for 
the previously underestimated values of elastic scattering 
time near the metal-insulator transition. 
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